JIB  nit 


O 


Department  of  the  Navy 
Scientific  Report 
Contract  N00014-76-C-0054 


HOMOGENIZATION  AND  LINEAR 
THERMOELASTICITY 


by 

Gilles  A.  Francfort 


SUDAM  Report  No.  81-4 


August  1981 


This  document  has  been  approved  for  public 
release  and  sale;  its  distribution  is 
unlimited. 


DIVISION 

OF 


APPLIED 

MECHANICS 


DEPARTMENT 

OF 


MECHANICAL 

ENGINEERING 


STANFORD 

UNIVERSITY 

STANFORD. 

UUFORNIA 

94305 


81  9  O'!  IGG 


UNCLASSIFIED _ 

SeCURITY  CLASSIFICATION  OF  THIS  PAGE  (Whan  Omt»  Bnlana) 


REPORT  DOCUMENTATION  PAGE 


.  REPORT  NUMBER 


title  Caid  Subtlllm) 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


T'S  catalog  number 


1  JIOMOGENIZATION  AND^INEAR  ^iERMOELASTICITY, 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

DIVISION  OP  APPLIED  MECHANICS 
STANFORD  UNIVERSITY 
STANFORD,  CALIFORNIA  94305 


II.  controlling  OFFICE  NAME  AND  ADDRESS  ^ - - 

OFFICE  OP  NAVAL  RESEARCH  (  /  { 

DEPARTMENT  OP  THE  NAVY  ''•U- 

ARLINGTON,  VIRGINIA  22217 


.  MONITORING  AGENCY  NAME  A  AOORESSflf  dlHannt  /ram  Controlling  Ollleo) 


10.  PROGRAM  ELEMENT,  PROJECT,  TASK 
AREA  A  WORK  UNIT  NUMBERS 

NR  064-525 


IS,  DISTRIBUTION  STATEMENT  (ol  Ihio  Roport) 


IS.  SECURITY  CLASS,  (ol  thit  nport) 


UNCLASSIFIED 


ISa,  DECLASSIFICATION/ downgrading 
SCHEDULE 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 


t7.  DISTRIBUTION  STATEMENT  (ot  th0  mbmtrmet  mitfd  tn  Btoek  20,  It  dm9rmt  from  Roport) 


18.  supplementary  notes 


19.  KEY  WORDS  (Contirtum  on  rovoroo  oido  it  noeomomy  ond  idonttfy  by  b'.ock  mmtbor) 


Dynamic  Thermoelasticity 

Composites 

Homogenization 


Periodic  Structures 
Asymptotic  Analysis 


JtPy^ABSTRACT  (Continum  an  rmrormm  »ldm  It  nocmoomry  and  tdontify  by  block  numbor) 

Me  study  homogenization  of  linear  dynamic  thermoelasticity  with  rapidly  vary¬ 
ing  coefficients,  using  a  semi-group  approach.  The  resulting  homogenized  problem 
exhibits  an  unusual  change  in  initial  temperature.  A  formal  asymptotic  analysis 
predicts  fast  time  oscillations  in  the  temperature  field.  These  oscillations  ex¬ 
plain  the  temperature  shift,  and  show  that,  for  our  problem,  weak  convergence  in 
time  is  the  best  convergence  that  one  czm  obtain.  / 


DD  1473 


EDITION  OF  I  NOV  AS  IS  OBSOLETE 


UNCLASSIFIED _ 

security  CLASSIFICATION  OF  THIS  PACE  (Whon  Doj^o  BntoroW) 


ABSTRACT 


We  study  homogenization  of  linear  dynamic  thermoelasticity  with 
rapidly  varying  coefficients,  using  a  semi-grov®)  c^proach.  The 
resulting  homogenized  problem  exhibits  an  unusual  change  in  initial 
ten^erature . 

A  formal  asymptotic  analysis  predicts  fast  time  oscillations 
in  the  temperature  field.  These  oscillations  explain  the  tenperature 
shift,  and  show  that,  for  our  problem,  weedc  convergence  in  time  is 
the  best  convergence  that  one  can  obtain. 

INTRODPCTION 

We  discuss  the  problem  of  "homogenizing"  the  equations  of  linecu: 
thermoelasticity  when  the  mechanical  and  thermal  properties  are  periodic 
and  rapidly  varying.  Following  Bensoussan  -  Lions  -  Papanicolaou  [1] 
and  Sanchez-Palencia  [7]  and  using  a  semi-group  approach,  we  show 
rigorously  that,  as  the  period  of  the  coefficients  goes  to  zero,  the 
solution  of  these  equations  converges  to  the  solution  of  a  related 
constant  coefficient  problem,  the  homogenized  problem.  Then  using 
a  formal  multiple-scales  method,  we  give  what  we  believe  to  be  a 
satisfying  interpretation  of  some  surprizing  features  of  the  results. 

Thermoelastic  behaviour  is  characterized  by  the  coupling  of 
hyperbolic  equations  of  motion  and  a  peuredxjlic  heat  equation.  This 
leads  to  several  interesting  phenomena  in  the  homogenization  process. 

Fast  oscillations  in  the  temperature  field  are  observed;  their  ampli¬ 
tude  remains  finite  as  the  period  goes  to  zero.  Thus  the  solutions 
C2U1  only  converge  in  a  weak  sense  in  time  to  a  slowly  varying  hcmogenized 


solution. 


Furthermore,  the  initial  data  for  the  homogenized  problem  are 
related  to  the  initial  data  of  the  inhomogeneous  problem  by  a  linecu: 
transformation  that  is  not  a  projection.  We  know  of  no  other  examples 
of  such  phenomenon. 

In  section  1,  we  formulate  and  prove  the  existence  of  a  homogenized 
thermoelastic  medium.  Section  2  contains  the  more  formal  eurguments 
euid  the  fast  oscillations  results,  which  are  at  the  root  of  the  ob¬ 
served  change  in  initial  data. 


1.  HOMOGENIZATION  OF  THE  THERMOELASTIC  PROBLEM 

To  cut  down  on  the  overwhelmingly  cumbersome  notations  that 
characterize  thermoelasticity,  we  will  place  ourselves  in  a  scctLar 
setting,  that  is  one  where  the  displacement  field  is  taken  to  be 
scalcur-vcilued.  Duvaut-Lions  [2]  show,  using  Korn's  theorem,  that 
this  is  no  loss  of  generality. 

We  consider  a  domain  Q  of  IR^  .  The  degree  of  smoothness  of 
will  depend  on  the  type  of  boundary  conditions  adopted.  We  will 
always  assume  that  is  smooth  enough  for  one  to  be  in  position  to 
apply  Rellich's  theorem  on  compact  imbeddings  of  SoboJev  spaces 
(Folland  [3] ,  Chapter  6) . 

We  will  refer  to  Y  =  II  ]0,  y°I  as  to  the  "reference  cell";  |y1 

i=l 

is  its  volume. 

If  E  is  a  smooth  hypersurface  dividing  Y  into  Y^  euid  ^2'  define 
aij(y),  6(y) f  P(y)  to  be  real  Y-periodic  functions, 

smooth  and  bounded  on  the  closure  of  Y^  amd  Y^  but  with  E  as  potential 
surface  of  discontinuity. 
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z 


Y 

Furthermore ,  ^ (y) ,  (y)  are  assumed 

to  be  symmetric,  strongly  elliptic  on  Y,  that 
is  that  there  exists  a  >  o  such  that: 

(1.1)  a^j(y)(resp.  X_(y))  j  ^ 

3(y)  and  P(y)  are  bounded  away  from  zero.  We  finally  choose  oi  such 
that  is  a  common  upper  bound  to  the  L^-norms  of  the  coefficients. 

We  extend  all  coefficients  to  all  of  by  periodicity.  Our  equations 
are  (Kupradze  [5]): 


(1.2) 


For  the  scUce  of  simplicity  we  will  only  ccnsider  Dirichlet  boundeury 
conditions  throughout. 

(1.3)  u^  =  0  ,  »  0  on  3R 


e 


3t" 


3 _  .X  3u- 

3x.  ^^ij  ^  3x, 
1  j 


a.  (g)  T  )) 


fi  A  51 
®  3t 


3 _  fX.  3t 

i  3 


®ij  l?lr 


And  for  initial  conditions,  we  will  have: 

(1.4)  u  (x,0)  »  f(x)  ,  (x,0)  =  g(x)  ,  T  (x,0)  *  lc(x) 


e  e 

Our  goal  is  to  study  the  behavior  of  u  and  X  as  £,  the  period, 
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goes  to  zero. 


We  define  H  to  be: 


(1.5)  H  -  H^(n)  X  L2(n)  X 


On  H,  we  define  the  operator  A^: 


\ 


(1.6)  Ag.  = 


1  3 _ ,  ,x.  3 _ . 


1  3 _ ,  ,x.„  ,x. .. 

P(£)  *ij  e  j 


1  ,x.^  ,x.  3 _ 

6.x.  ®ij^e^  j^e^3xi 
'e' 


1 


6,5, 


with  domain 

(1.7)  D(A_)  =  {u  =  (U,  u.  ,  X)GH^(n)  X  L,(n)  X  (fl) 

c  t  o  2  .  o 


such  that  A_o  G  H  in  a  weak  sense  } 


Then  the  following  proposition  holds: 


Proposition  1.1 


Ag.  generates  in  H 


a  strongly  continuous  semi-group  of  operators 


Sg.(t)  such  that  : 

(1.8)  I  I  Sg(t)  1  I  1  (V  t  >  0  ) 
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Proof 


We  first  consider  for  a  fixed  £  the  norm 


(1.9) 


where  ~  denotes  con^lex  conjugate. 

In  view  of  the  properties  of  the  coefficients,  |  |  g.  is  indeed 
a  norm  on  H,  equivalent  to  the  natural  Sobolev  norm  on  H,  noted 
II  II  ,  that  is,  if  U  is  in  H: 


(1.10)  Qt  II  u||^  £  |u|e  l“*^ll  u||^ 


In  the  norm  |  |  g  »  generates  a  semi-group  of  contractions. 

00 

Indeed,  the  domain  D(Ag.)  is  dense,  since,  though  (J2)  functions  do 

00 

not  belong  to  it,  fimctions  whose  normal  derivatives  cure  0  on 

the  only  possible  surfaces  of  discontinuity  for  the  coefficients  (i.e. 
the  ^-scaled  versions  of  I  in  each  of  the  cells  mcUcing  up  n)  do  belong 
to  the  domain  D(Ag.)  .  Checking  that  Ag.  is  closed,  that  the  range  of 
(1-Ag)  is  H  itself  and  that  Ag.  is  dissipative  offers  no  special  dif¬ 
ficulties  (see  Frauicfort  [4]  for  full  details).  Note  that  the  measure 
of  the  dissipation, 

(1.11)  Re(A£U,  U)  -  -  Re  (/rj  1  ’  “l^^lL(JJ) 

3  i  2 

(in  view  of  the  properties  of  the  's),  is  precisely  the  physical 
dissipation  due  to  heat  fluxes  through  the  boundary. 


13 


The  result  then  follows  from  the  application  of  Lumer-Phillips's 
theorem  (Yosida  [8],  Chapter  9).  Therefore, 

(1.12)  lSe{t)ule_<  |ulg.  for  any  U  in  H 

and  thus,  using  (1.10), 

(1.13)  II  Sg.(t)u|l  la’^liuil 


which  conpletes  the  proof. 

We  now  leave  the  time  dependent  formulation  and  examine  the 
behavior  of  the  resolvent  of  (Ag.)  cis  £  goes  to  0.  At  the  end  of 

this  section  we  will  reintroduce  the  time  dependence  by  using  some  basic 
properties  of  semi-groups. 

It  is  a  direct  consequence  of  (1.8)  (Yosida  [8],  Chapter  9)  that 
the  right  half  coiqplex  plane  belongs  to  the  resolvent  set  of  A^,  for 
every  £.  Let  us  consider  F  =  (f,  g,  k)  to  be  an  element  of  H.  Taking 
X  to  be  real  strictly  positive,  we  have  the  following  string  of 
eq\ii  valences : 


(1.14)  (A£.)F  =  Ug.  ,  (U£  =  (u 


e  e 
u*.  ,1  )) 


Xu  "  ■  f 


P(f)g 


Leif) 


,x,  9t 


3u^ 

Xn  t 


&(4)k 


i  3  3. 


6 


o 


(1.15) 


£  £ 
Xu  - 


I  X^P(i)u^  -  -Ij-  (a_{^)  (|^  -  ujXjxE))  =.  p(X)  + 


ij  e  '■55r  j'£'  "  -£' 


g) 


.X  e(|)T<^  -  |_(X^^,|)^).  X.^^(|,a,,|,,^.  6(|)k  *  .ij  (1)0*3  (fill- 


C  £ 

The  last  two  equations  (1.15)  have  a  unique  solution  v  =  Xu  ,  T 
1  2 

in  (H  (Q) )  ,  since  the  Dirichlet  form  dp  defined  as: 

o  ^ 


(1.16)  d,«v^T^,(;,m)  =  i  4  |S-  dx  . 


X/„  P(f)v^?  d*  -  /jj.a.  .(Ijajdlt*^  li;  ^  *  ^S1 


*  X/jj  8(|)t^n  d*  .  4  "  *■ 


1  2 

is  strictly  coercive  on  (H  (J2) )  ,  in  view  of  the  properties  of  the  coef- 

o 

ficients. 

£  £ 

If  we  manage  to  find  a  limit  for  u  ,  T  as  £  goes  to  zero,  then 
going  back  up  through  the  string  (1.14)  will  enable  us  to  get  the 
limit  of  R^(A£)F. 

Performing  the  limiting  pi;ocess  in  (1.14)  is  the  task  of  the 
homogenization  method.  Rather  than  going  through  all  the  lengthy  details 
of  the  argument,  we  merely  mention  the  different  steps  that  were  per¬ 
formed,  underlining  only  the  ones  that  are  not  totally  stemdcurd.  For 
further  detedls  the  reader  is  to  refer  to  Bensoussan-Lions-Papanicolaou 
[1] ,  Chapter  1*  especially  Sections  3,  9,  and  13,  or,  for  our  problem  to 
Francfort  [4] . 
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First  one  shows  that  Ug-  and  are  bounded  in  which 

immediately  implies  the  existence  of  a  wecdcly  convergent  subsequence 
1  2 

in  (H^(n))  converging  to  (UjT) .  Since  we  ultimately  show  that  emy 
convergent  subsequence  converges  to  the  same  limit  we  do  not  distinguish 
between  the  sequence  and  siibsequences  of  this  sequence. 

Then,  defining 


^CJi  =  - stress, 

X]  <JX^ 


(1.17) 


1  ij  ^  oXj 

1]  £  j  ^  9xj 


the  heat  flux. 


it  is  ecisy  to  conclude  that  these  quantities  converge  weakly  in 
to  which  in  turn  satisfy: 

C-  2 

pX^u  -  - —  =  P(Xf  +  g) 

1  3-i 

(1.18)  < 

9<.  _  ^4= 

HKt:  -  ^  +  Xv  =  3lt  +  a.  .  a.  ^  , 

L  ^  ID  D  9x^ 

where  ,  from  now  on,  will  denote  the  Y-average  dy. 

It  remains  to  determine  and  v.  This  is  the  core  of  homogen i 

zation.  TO  this  effect  we  define  Xj^  ^Y)  »  ^^® 

periodic  solutions,  up  to  a  constant,  in  H^(Y)  of: 


8 


(1.19) 


3 


a 

^  (kij  (y)  ^) 


a  3'!' 

•3?"  '*13''”  ■55"’ 


= '  -577 


H'  can  be  considered  as  non  standard  with  respect  to  the  "classical" 


caise .  The  functions : 


(1.20)  «£  =  -  e  ^  ^  ®k^£^  ' 


satisfy: 


(1.21) 


3we  a** 

.  X.  _ k  3u) 

•'n  *ij^e^  3Xj  3xj^ 

..  9z,^ 


dx  =  0 


/o  X.  .  (f)  ^  ^  dx  *  0,  for  any  ^,\i  in 
•'ii  1]  E  3xj  ox^ 


Taking  u  and  U  to  be  C^Ca)  functions  and  making  use  of  (1.16),  (1.21) 


we  have: 


(1.22) 


X,  ^k  3  (ii>u  ) 


d_((Xu  ,  T  ),(W«k'UZjc^)  ^ij^F  ^ 


3z^  £ 

X.  .(|)P^  dx  -  /  P(f)  (Xf  +  g)  ^  W^  dx 

11  t  3xj  3x^  “  ^ 


.;n(B(|)k  1  a,.(|)a,(f)||-)u-z^dx 
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In  (1.22),  we  have  in  essence  siibtracted  from  the  variational 
formulation  of  (1.15)  appropriate  expressions  equal  to  0  in  order  to 
eliminate  products  of  weak  convergences. 

It  is  then  possible  to  go  to  the  limit  in  (1.22)  in  a  way  identical 
to  Bensoussan-Lions-Papanicolaou  [1],  Chapter  1,  Section  3.  Upon 
performing  this  limiting  process  and  come  out  to  be: 


^i  .  9u  _  ,  _  a.  ^ 

^ij  “  ^kj  ^  ■  ^ij  “j  "  j  9yj^ 


(1.23) 


=  a  ,  .  9t 

“^i  “  ij  “  \j  3x^  9x. 


Determining  V  requires  some  extra  effort  and  the  use  of  f.  One 
e 

defines  g  to  be: 


(1.24)  g^  =  1  +  e  4'(|) 


then  it  satisfies,  for  ciny  in  (0)  : 


,  X  3g  3iij  ,  ,x,„  ,x,  3ai  , 

(1.25)  3x.  -9^- 


Going  through  the  same  procedure  as  in  (1.22)  but  with  U  equal  to  0 
£  £ 

auld  Wj^  replaced  by  g  ,  we  determine  v  to  be: 


af  ,  3u  ^  ,  9'i'  ^  T 

(1.26)  V  =  a.  a.  -  a  ^  t —  +  (a  a,.- —  )  i 

13  3  3y  3x^  k]  33y^ 


10 


o 


Defining  ,  A^,  Y^,  O  to  he 


^  -  "kj 


1]  13 


3Xi 
k 


an? 


(1.27)  / 


B.  “a..oi.-a..„ 

1  11  1  11  3yj 


30. 

X .  .  =  X .  .  -  X  -~ 
11  11  k]  3y 


-  a.  .ot,  -  A. 

1  11  1  1 


it  can  be  shown,  using  (1.19),  that  a^^^  and  X^^  are  symmetric  positive 


definite,  hence  invertible,  that  A.  emd  B.  are  equal  and  that  O  is 
positive. 

We  set: 


(1.28)  “  ^ik  ^  ■  ^ik  ®k 


Recalling  (1.18),  (1.23),  (1.26)-(1.28)  yields: 


(1.29)  < 


‘1^  -  “j  i-  ’  ‘ 


(^  +  o) Xt  -*  X.  ,  ~ — -  +  Xa,  .ct .  -K— —  =  hk  +  i (y) *1  •  (y) 5  I 

11  ij  1  11  1 


1  1 


11 


and,  in  view  of  the  properties  of  the  a^j's  emd  the  Dirichlet 

1  2 

form  cissociated  to  (1.29)  is  strictly  coercive  on  (H  (H) )  ,  hence 

1  2 

(1.29)  admits  a  unique  solution  in  (H  (H) )  .  Then,  using  (1.14), 

o 

we  end  up  with  the  following  proposition: 

Proposition  1.2 


R^{A^)F  converges  weakly  in  (H^((2))^  to  the  unique  solution  in 
(H^(fl))  ^  of; 


Xu  -  u^  =  f 

(1.30)  ixpu^  --ij  <1^ 


8  T  t 

X(-5  +  0)T  -  X^  +  a  c.  ^ 

i  i  1 


-  ^  *  ''i  % 


We  then  define  A  to  be: 


(1.31)  A  = 


1  ,  3 

-  ij  3x. 9x. 
P  IT 


$  +  a 


i.  .  ot  — 
1] 


-  i  a.  a.  I — 
^  1]  3  9x^ 


I  .  o  i33x^3Xj 


It  is  simply  a  matter  of  reproducing  the  proof  of  proposition  1.1, 
but  with  constant  coefficient  this  time  to  show  that  A  generates 
a  semi-group  of  operators  S(t)  such  that: 


12 


(1.32) 


II  s(t)  II  <_  a' 


for  euiy  t  >  0 


-1  •  '-1  .  . 
Renaming'  ^  fhe  maximum  of  ^  and  ct  ,  we  deduce  from  proposi'tron 

1.2  euid  (1.32)  the  following  corollary. 

Corollary  (1.2) 

R^  (Ag.)  F  converges  wecJcly  in  (H^  (f2) )  ^  to  R^^  ( A)  F  where ; 


(1.33) 


+  Y.  1^ 

X  3x. 

(f»  g»  — = - -  ) 

3  +  a 


Now,  (1.8)  implies  that,  for  any  U,  there  is  a  bounded  subsequence 
of  Sg.(t)C  that  converges  we2ik-*  in  ^  element  of 

1^0(184.  ,  H) .  This  is  a  direct  consequence  of  the  sepeurability  of 
Lj^OR+ ,  H)  and  of  Banach-Alaoglu's  theorem  (Rudin  [8],  Chapter  3). 

Still  identifying  a  sequence  with  its  sxibsequences ,  we  get  that,  for 
any  V  in  H, 

00  «  QO  . 

(1.34)  e“'^^{Sg.(t)U,  V)  jj  dt  ^^®(t)  ' 


where  (  f  )jj  Is  the  natural  inner  product  on  H. 

But  the  resolvent  of  the  generator  of  a  semi-group  applied  on  a  vector 
U  is  equal  to  the  Laplace-transform  of  the  semi-group  acting  on  U 
(Yosida  [8],  Chapter  9)  thus; 


(1.35) 


e”^^(Se(t)U,  V)^dt 


(Rj^(A£)U,  V)jj 
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•*«* 


which  itself  converges  to: 

(1.36)  (R;^(A)y,  V)jj  =  e"^^(S{t)y,  V)jjdt 


Since  V  is  arbitrary,  we  finally  get,  using  the  uniqueness  of  Laplace 
transforms  of  scalar  function  that: 


(1.37)  G(t)  “  S(t)y  (t  >  0) 


We  have  proved  in  this  section  the  following  theorem: 

Theorem 

The  generalised  solution  of  (1.2)  with  Dirichlet  boundaury  conditions 
and  initial  conditions  (f ,  g,  k)  in  H  converges  weak-*  in  ,  H) 

to  the  generalized  solution  of  : 


(1.38)  < 


'■pii 


a.  .  ( 


^2 
d  u 


ij ' 3x. 3x . 
1  3 


-  01  , 


V, 


ij  3x .  3x . 
13 


3T^  , 

X 


-  a. .  a . 
13  3 


a2 

0  u 
3t3x. 


with  Dirichlet  bound6u:y  conditions  and  initial  conditions 


(1.39) 


(f,  g,  — -  ) 

?  +  a 


14 


Before  concluding  this  section,  let  us  emphasize  once  more  the 
rather  unusual  change  in  initial  temperature  in  (1.39). 

2.  FAST  OSCILLATIOtlS  OF  THE  TEMPERATURE  FIELD 

Since,  through  a  weak-*  type  of  convergence  a  rapidly  oscillating 

function  (like  e^^^  )  goes  to  0,  it  is  fairly  natural  to  expect  a  ^ 

e  .  _e 

dependence  of  u  and  ^  ,  This  kind  of  problem  is  most  easily  addressed 
using  asymptotic  expansion  techniques.  We  have  already  mentioned  the 
semi-heuristic  chcuracter  of  this  section  so  that  we  will  not  dwell 
on  the  restrictions  to  the  problem  that  would  make  the  aurgrsnent  totally 
rigorous . 

£  £ 

Recalling  (1.2)  we  now  suppose  that  u  and  T  are  functions  of 
t  and  6»  t/E;  3^  becomes  ^  We  then  Laplace  transform 

(1.2)  with  respect  to  both  t  and  5,  the  dual  variables  being  respectively 
^  and  u  •  From  now  on: 

will  denote  the  t  -  Laplace  -  transform 
will  denote  the  6  -  Laplace  -  tremsform 
''  will  denote  ~  or  ^ 

In  order  to  be  able  to  perform  these  transformations  we  need  to 
impose  initial  conditions  on  both  t  and  6.  we  will  set: 

u^(x;  0,  a  f(x)  ,  u^(x;  t,  0)  *  p(x,  t) 

£  -  e 

(2.1)  3^^**  0,  6  )  a  g(x)  ,  t,  0)  »  q(x,  t) 

T^(x;  0,  6)  »  k(x)  ,  (x;  t,  0)  ■  0(x,  t) , 
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where  £,  g,  k  are  as  before  emd  p,  q,  6  are  unknown.  We  get: 

r 

^  f  +  I  -  ip)  + 


(2.2)  / 


1  (U2u^  UP  -  ^)  }  =  Ij-  (a.^(|)  (  f  -  -  a.  (|,t" 


)) 


e(|){(;x^  -  ^)  -  |(ux"  -  ®>>  =  It  ’ 

.»  i  j 


.X.  ,x.  1-3  •r'^e  fv  1  9  >'£  /V.  1 

-  a  (£)a  (£){^  (iu  -  -)  +  -  ^  (pu  -  p)} 

i  X 


£  e 

We  seek  an  eaqpansion  of  u  aind  T  in  the  form 


(2,3 


i 

Z  £  'i^(x,  y,  t,  6 ) 


=  Z  £^  T^(x,  y,  t,  6 )  ,  where  y  =  ^  . 

The  dependence  of  the  ' s  and  T^' s  on  y  is  taken  to  be  Y-periodic. 
This  is  always  what  is  assumed  when  performing  double  scaling  in  space 
in  problems  related  to  homogenization. 

We  also  need  to  control  the  faist  time  behaviour  of  u^  and 
Since  we  would  like  them  to  be  oscillating  in  <5,  or,  at  leeist,  to  be 
such  that 


(2.4)  ““  T  -►  +  «>  ^  /  u^(x,  t,  y,  6)  dfi  (respectively  T^) 

o 

exist  and  be  finite,  we  are  led,  through  Wiener's  Tauberian  theorem 
(Rxidin  16],  Chapter  9)  to  suppose  that: 
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(2.5)  lim  yu-  (respectively  exists  and  is  finite, 


wd  we  will  furthenoore  assxxne  that  this  linit  is  to  be  taken  point- 


wise  in  X  and  weakly  in  H‘‘’(Y)  with  regard  to  the  y  dependence. 

With  these  considerations  in  mind  we  can  proceed  to  replace 
3  1  3  £  e 

by  and  u  and  T  by  their  expansions  in  (2.2) 


3x. 


-2 

We  obtain  two  "series"  in  ciscending  powers  of  £  starting  at  £  ; 

we  identify  the  factors  of  each  of  these  powers  to  0  ,  one  after  another. 
_2 

As  factor  of  £  we  get: 


3^S 


(2.6)  / 


p(y)  (u  Ug  -  UP  -  q)  »  ^  ^  ^ 


3T 


-  p)  =  0 


Since  the  Dirichlet  form  cissociated  to  the  operator 


(2.7)  D  =  P(y)U^*  -  1^  (a..(y)  ) 


is  strictly  coercive  on  the  subspace  of  H''’(Y)  consisting  of  Y-periodic 

functions,  the  first  equation  of  (2.6)  has  a  unique  solution  in  H  (Y) ; 

therefore  ^  solution.  Hence 

P  y2 


(2.8)  .  a 

^  V 

But,  in  view  of  (2.5) ,  (2.8)  implies  that  q  ■  0  thus  u^  is  equal  to 
t  and  does  not  depend  on  y.  Inverting  (2.8)  we  get  that  u^  does  not 
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o 


depend  on  6  either; 

(2.9)  Uq(x,  t)  -  p(x,  t) 


Then  from  the  second  equation  of  (2.6), 


V  V 


(2.10)  T  -  T  (X,  M)  , 
o  o 


since  the  only  periodic  solution  of  that  equation  is  a  constant  with 
respect  to  y. 

As  factor  of  £  ^  we  get,  using  (2.8),  (2.10): 


r 


(2.11) 


Du,  = 


Si(y]P^o  3 

3yi  '  ’3y“  '^ij 


(a.  Jy)a.(y))  c 


6(y)  (UT^  -  0)  -  3^  >  +  37^  ^ 


3u^ 

-Ua^.(y)a^(y)  ^ 


Defining  to  be  the  unique  periodic  solutions  in  H^(Y)  of 

V 


(2.12)  / 


2  u  3  ^i 


3yv 


(y) 


we  obteiin  from  the  first  equation  of  (2.11) 


(2.13)  u. 


3u  Li 

u  o 

-Xj  ^ 


Then,  integrating  the  second  equation  of  (2.11)  with  respect  to  y  and 
defining  and  to  be  the  amalogues  of  Y^  and  o  for  and  as  in 
(1.27) , 
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■gg  +  Y_U  |i_  (def.) 


1  ax^ 


where  denotes  the  Y-average  dy 

We  introduce  and  to  be  the  unique  periodic  solutions  in  H^(Y), 
k 

up  to  a  constant,  of: 


(2.15) 


- <’'ij  <y>  3^> 

■157  ' 


■K  k 


■  -  <niW“jW  357  -  ^ 
-  (a..(y)“.(y) 

i 


The  equations  (2.15)  are  well-posed  since  the  Y-average  of  the  right- 
hamd  members  do  vanish  by  definition  of  euid  0^  . 

Recalling  0^ (y) ,  we  obtain  for 

(0  X  =  -  —  0  fv)  +  H*^r)  +  an  arbitrary  function  of  x  onlj 

1  M  j  ^  3x.  j  ax.  U 


Finally,  as  factor  of  e  ,  we  get; 


^  ^  vu  - 

P(y)  ^  -  y  +  2C)iu^)  +  DU2  » 

2  2 

as  as,  a 

+  a.  .  (v)  ^  +  a  (v)  -  -  a.  (y)*^  •  (y)  .r  *  -  j  (e.  .(y)oi.(y)"^,) 

^  ax^ax.  ij'^^ax  ay  r^'ax^  ^  ij  :  ^ 


o 


Wh 


6(y)  (5^„  -  ^  ^  (y)^  )  *  |j-(X^ .  (y)^  ) 


(2.17) 


3^T 


-  a 


ij(y)0tj(y)  1^(52^)  -  ^j(y)“j(y)^(Uiii)  -  a.,(y)a.(y)^(uu2) 


Vte  integrate  both  equations  of  (2.17)  with  respect  to  y;  making 
use  of  all  the  previous  results  of  this  section,  we  get: 


(2.18) 


^P(?^p  -  5f  -  g)  +  ^  .1^ 


2  i  j  3x .  3x . 

i  D 


a'!,  a^! 
i-D  3 


!(.?  .  k)  .  y^  6^,  .  .  1^ 

J 


u  x..(y>^  ^ 


•  •(y)3r^  +  vix.  .(y)|S_^ 

1]  ^  3yj  3x^3x^  13  ^  3yj  3x^ 


13  3  3«i  i^(y)  x,y,  3»i 


^  ,.2_  _ ^*^3 


y  a^j  (y)oij  (y)-g5J-  '  U  (y)aj  (y)-g^  » 


3x. 

1 


where  a^^^,  are  to  Xj^  ^Jid  4^  what  a^^  and  are  to  Xj  and  f  in  (1.27)  . 

We  now  consider  the  limit  of  (2.18)  as  y  goes  to  0.  The  following 
result  holds: 

Proposition  2.1 

Xk^f  4*^,  y  yH^  go  respectively  to  Xj^»  'J*  »  0  and  0  strongly 
in  H^(Y)/K  as  y  goes  to  0.  Hence  a  ^,ci  ^,Y  ^,0^  go  to  a  ,cx  ,Y.  ,  <7. 

ij  J  1  ij  J  1 
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The  proof  of  this  proposition,  which  involves  some  basic  esti¬ 


mates  in  H^(Y)/3R  will  not  be  given  here;  refer  to  Francfort  [4] 
for  the  details. 

Proposition  2.1  together  with  (2.5)  is  exactly  what  we  need 
to  perform  the  limiting  process.  Upon  doing  so,  we  come  up  with  a 

A 

set  of  two  equations  for  p  and  H  which,  together  with  the  limit  of 
(2.14),  can  be  interpreted  in  the  time  dependent  domain,  p  and  0 


satisfy: 


(2.19) 


0!x,  t) 


(3  +  cr)n(x,  t)  -  (x,  t) 

*i 


1 


T  3 


3^ 


■  3n 

3^2  ~  ^ij  3x^3Xj  “  *ij  j 


-  a.  .a. 


3_E_ 


ij  j  3t3x^ 


3f 

-  "Sc  +  Yi  3x^ 

p(x,  o)  =  f,  (x,  o)  *  g,  r|(x,  o)  =  - 

B  +  a 


where  U  is  the  limit  of  IT  as  U  goes  to  0. 

It  is  clear  that  ri(x,  t)  can  be  identified  with  T{x,  t)  ,  the  homo 

genized  temperature  field,  and  p(x,  t)  with  u(x,  t) ,  the  homogenized 

displacement  field.  Replacing  0  by  its  value  in  (2.14)  we  also  obtedn 

e 

w  e^qjression  for  the  leading  term  of  the  eisyirptotic  expantion  of  i  , 
that  is  its  <S-Laplace  transform  satisfies: 


(2.20) 


X  .i  2i!.  ^ 
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a  . 


(2.20)  is  not  invertible  in  general,  but  the  following  proposition  holds 
Proposition  2.2 

emd  go  to  zero  eis  u  goes  to  + 

The  proof  of  this  last  proposition  uses  the  same  estimates  eis 
the  ones  that  establish  proposition  2.1. 

Propositions  2.1  and  2.2  enable  us  to  conclude  that,  as  u 
goes  to  0,  goes  to  n,  whereas  as  U  goes  to  +  »,  goes  to  0. 

In  a  time  dependent  context,  these  facts  translate  into  statements 
on  the  behavior  of  T  near  infinity  cmd  near  the  origin. 


(2-21) 


lim  j  _ 
0  -+■  +  00 


^lim 


6  -»•  o  + 


T:{x,t) 

0(x,  t) 


The  second  equation  of  (2.21)  is  consistent  with  our  self  imposed 
6  -  initial  conditions.  The  first  equation  shows  that  the  fast 
oscillations  of  the  leading  term  of  the  asymptotic  expansion  of 
are  centered  about  T(x,t) ,  the  solution  of  the  homogenized  problem. 


CONCLUSION 

Numerical  evidence  corroborates  the  results  of  section  2  and 
confirms  that  feist  oscillations  are  indeed  the  phenomenon  leading 
to  this  unusual  change  in  initial  data  [4] . 

If  seeking  a  more  physical  explanation,  one  could  examine  the 
entropy  associated  with  the  problem; 
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8u^ 

3x. 


1 


It  is  fairly  straightforward  ,  using  the  results  of  section  2  and 

some  of  the  steps  performed  there,  to  show  that  there  is  no  fast 

dependence  in  time  of  the  space  average  of  the  leading  term  in  the 
e 

e^^ansion  of  s  ■  That  the  macroscopic  entropy  of  the  body  is  a 
slowly  varying  quantity  appears  to  be  a  sound  idea  emd  does  fit 
our  physical  intuition.  A  fast  oscillation  in  the  temperature  field 
is  the  effect  that  balances  the  space  oscillations  of  the  strains 
due  to  the  inhomogenities  of  the  coefficients  euid  allows  the  entropy 
to  evolve  slowly  at  its  own  pace.  In  this  respect  the  unusual 
initial  change  in  temperature  is  needed  to  insure  that  no  fast 
change  in  entropy  is  taking  place  at  time  zero. 

In  contrast  with  other  fast  oscillation  type  problems,  the  "phase" 
of  the  oscillations  in  not  arbitraury  but  perfectly  determined.  It 
also  appears  that  a  geometrical  optics  type  ansatz  in  place  of 
(2.1)  will  fauil  since, if  the  solutions  of  (2.12)  are  sums  of  terms 
of  more  than  one  frequency  in  <S,  the  fast  oscillations  need  not  be  periodic 
in  6. 

To  conclude  this  study  let  us  point  out  that  choosing  the 
entropy  as  the  nat\iral  vauriable  in  place  of  the  tenperatxire  introduces 
space  derivatives  of  the  third  order  and  thereby  prohibits  a  rigorous 
analysis  of  the  type  performed  in  section  1.  A  perturbation  analysis 
using  double  scaling  is  feasible  but  eventually  leads  to  reintro¬ 
ducing  the  temperature  field  as  the  proper  variable. 
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